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Abstract 

A  closed  expression  is  given  for  the  distribution  functions  in  the 
classical  theory  of  statistical  equilibrium.  This  expression,  when  suitably- 
expanded,  permits  one  to  obtain  rather  simply  the  general  terra  in  the  virial 
expansion  of  the  distribution  functions  and  of  the  equation  of  state. 
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1.  Introduction 

The  object  of  the  present  note  is  to  give  an  explicit  solution  of  a  system 
of  equations  introduced  by  Bogolubov'-  -■  for  the  distribution  functions 
in  the  classical  theory  of  statistical  equilibrivm.  This  solution,  when 
suitably  expanded,  permits  one  to  obtain  rather  simply  the  general  term  in 
the  virial  expansion  of  the  distribution  functions  and  of  the  equation  of  state. 
Although  our  solution  is  given  in  closed  form,  it  must  be  considered  as  being 
of  formal  nature,  since  it  involves  operations  of  differentiation  and  of 
substitution  of  variables  to  be  performed  on  functionals. 

We  consider  a  system  of  N  identical  monatomic  molecules  enclosed  in  a 

macroscopic  volume  V.  It  will  be  assimied  that  both  N  and  V  tend  to  infinity, 
in  such  a  way  that  their  ratio  remains  constant 

(1)  I  =  c- 

The  symmetric  distribution  functions  F  (x^,...x  ),  normalized  so  that 


(2) 


Jio^/..-i,/='^^'-'^^'^^-^=° 


V->«,  .   .y    -y 


satisfy  the  infinite  syston  of  covrpled  eqiiations 
SF  au  r    S$(|x-x       I) 


s  =  1,2...; 
a  =  1,2,5. 


1  1 

where  J  is  the  intermolecular  potential,  U^  =•  g  ZI  1(  \'\'\  1^'  ^  ^  kf  * 

The  infinite  system  of  equations  (5)  has  been  derived  by  Bogolubov 

from  the  customary  canonical  distribution  of  Gibbs.  An  alternative  derivation 

presents  some  interest,  since  it  may  provide  an  independent  foimdation  for 

equilibrium  statistics.  One  first  obtains,  from  the  Liouville  equation,  the 

well-known  hierarchy  of  infinitely  many  equations  for  the  time-dependent 
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distribution  functions '-^-' *'--'* '-1*^-'.  One  then  seeks  time-independent  solutions  which 
are  Maxwellian  in  their  velocity  dependence:  for  their  cooi^iinate  dependence 
the  system  (5)  obtains. 

The  infinite  system  (5)  is  clearly  equivalent  to  a  single  functional 
differential  equation  for  the  generating  functional 


m  =r  7r/.../vx,,...,x^)'X(x,)...x(x3)dx^.., 

from  which  the  distribution  functions  can  be  obtained  by  functional 
differentiation 

(5)  ^s^\'""-s^=w^hmrj 

^  '   X  =  o 

The  functional  F  is  defined  for  test  functions  X  which  tend  to  zero  sufficiently 


fast  at  infinity.  The  equation  satisfied  by  F  is 

3    6F       r  ^i^Hx-yl)    52p     ,^  ,  ,   , 
(6)     —-  ^rJf^s   +   M    —7 rV/^\r-V/„\  (X  (y)+c)dy  =  0  . 


The  physically  meaningful  solution  of  (5)  or  (6)  is  known  in  the  limit 

** 
c  ->  0,  of  zero  density  : 

-txU 
(7)      F°  =  1,  F°  =e   ^     s   =   2,5,...   . 

The  equations  can  be  solved  in  power  series  of  c,  starting  from  (?)  as  a 
zero-order  solution.  The  first  few  terras  have  been  calculated  and  shown  to 

^ _ 

Equations  of  this  type  have  been  considered,  e.g.,  by  Bazarov '•^-' . 

** 

Here  and  in  the  following  we  shall  simply  denote  by  a  superscript  zero 

quantities  referring  to  the  zero  density  limit.  They  will  be  related  to 

each  other  by  equations  similar  to  those  relating  the  corresponding 

quantities  without  superscript.  Thus  one  has,  together  with  equations 

(10)  and  (12)  below,  the  equations 


,0  _    e^O 

F^  =  e"'  —1, 


,0   ^G""  ,      gc   6G"' 


Quantities  with  a  superscript  zero  are  considered  known. 
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reproduce  the  well-known  virial  expansion  .  However,  the  form  of  equation  (6) 
suggests  that  one  may  obtain  a  complete  solution  by  taking  as  independent 
variable  the  function 

(8)  X  (y)  =  X  (y)  +  c 

since  the  equation  then  reduces  to  its  form  for  c  =  0.  This  is  not  quite 

correct,  of  course,  because  X  is  not  an  admissible  test  function  for  F 

(it  does  not  vanish  at  infinity) ^  some  transformations  inspired  by  the  theory 

of  Ursell  function  must  be  performed  first. 

2.  Formal  expression  for  the  distribution  functions 

To  canry  out  the  transformations  mentioned  in  the  preceding  section,  we 
introduce  a  new  functional 

oo  ,   / 

(9)  G[X]  =11    JT  \   Gs^V---^)  X(xi)...X(x3)dx^...dx^ 

8=1     J 

related  to  f[X]  by 

(10)  F  =  e^  -1   . 

The  functions  G  generated  by  G  are  related  to  the  distribution  functions 
F  in  the  same  way  that  the  customary  Ursell  functions  are  related  to  the 
Boltzmann  factors 

Fl(x^)  =  G^(x^) 

(11) 

F^{x^,x^)   =  G^{x^,x^)   +  G^{x^)G^{x^) 

etc. 

See  the  papers  quoted  in  [ij  and  [5]. 

The  connection  between  eq\iations  like  (lO)  and  (ll)  has  been  pointed  out 
before  by  Mazo  and  Zemach'-  -^ .  These  authatts  also  show  how  (lO)  can  be  used 
to  p3X)ve  quite  siniply  the  property  of  the  G  functions  of  vanishing  for 
lai^e  sepai 
separation. 


k  . 


The  equation  (6)  can  be  quite  easily  transformed  into  an  equation  for  G,  or 
rather,  for 

00 


(12) 


^b'^  -wh'^*i:^hj  =.,,(x,x^,...,x^)X(v...\(xjd^...ax^ 


It  takes  the  form 

(15)   ^  h[x,  X]  *  .  /  ^  (X  (y)«)(5i|^  -  H[x,t]H[y,X])dy  =  0  . 

The  advantage  of  this  transformation  is  that  X  +  c  is  an  admissible  test 
function  for  H  x,  Xjj  since  the  functions  generated  by  H  are  of  the  Ursell 
type,  and  vanish  themselves  when  the  coordinates  are  far  away  from  x,   so  that 
the  test  function  does  not  have  to  vanish.  The  substitution  (8)  is  now 
permissible.  More  generally  one  can  easily  verify  that 
(li.)    h[x,X]  =aH°[x,a06+c)] 

is  a  solution,  a  being  an  arbitrary  constant. 

The  constant  a  cannot  be  chosen  eqixal  to  unity,  since  the  normalization 
condition  (2)  implies  that 

(15)  G^(x)[  =  F^(x)]  =  1 

and  therefore,  from  (12)  and  (l^),  that 

(16)  aH°[x,ac]  =  1. 

For  the  purpose  of  identification  with  known  resxilts  one  ceui  show  quite 
easily  that  the  condition  (l6)  for  a  can  be  written  as 

(17)  af:  nb^(ac)°-^  =  1 

n=l 

where 
(18) 


^1  =  ^'  \=^|^n(W-^n)^S 
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are  the  well-known  cluster  integrals  L-^-' .  Therefore  the  product 

(19)  ac  =  z 

is  the  activity  and  (I6),    or 

(20)  c  =  zH°[x,z]    , 

expresses  the  relation  between  density  and  activity.  Equation  (l4)  can  now  be 
written 


(21)      h[x,X]  =f  H°[x,|X  +  z]  . 


Equations  (9),  (12),  (20)  and  (21)  together  describe  the  desired  solution. 
The  solution  appears  expressed  in  terms  of  the  activity  z;  if  one  desires  an 
expression  in  terms  of  the  density,  one  must  first  invert  eqviation  (20). 
5.  Expansion  to  first  order  in  the  density 

As  an  example,  we  shall  now  use  the  obtained  solution  to  calciilate  the 
distribution  function  F  neglecting  terms  of  second  order  in  the  density. 
From  (20) 

(22)      f  =  H°[x,z]  =  1  +  z 

Equation  (22)  shows  that 

(25)      z  =  c  +0(c^) 

so  that,  in  the  following,  we  will  have  to  retain  only  terms  of  first  order 


f  G°(x,x^)dx^ 


+  0(z^). 


in  z.   Clearly,  from  (21), 


(24)      G2(x^,X2)  =  -5^ 


^^r^] 


zx2 ""  ri'^J 


6E" 


(|) 


-m 


and 
(25) 


t-^ 


G°ix^,x^)   +  z  I  G°(x^,X2,x^)dx^  +0(: 


-)    . 
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Finally,  substituting  in  (2ij-)  the  results  (22)  and  (25), 
(26)    G^Cx^^Xg)  =  G°(: 


^l'^2^  {-^'^   J  G°(x^,x^)<i)c^-c  f  G°(x2,x^)dx^^ 


j'>i'-2' 


x^)dx^+  O(c^). 


The  function  F  is  given  by 

(27)     F^Cx^^x^)  =  G^Cx^^Xg)  +  1 

It  is  quite  easy  to  see  that  the  above  approximate  result  agrees  with 
previotisly  given  expressions.  As  it  is  shovn,  for  instance,  in  Reference  1, 
this  result  can  be  used  to  obtain  the  virial  expajision  of  the  equation  of  state 
to  second  order  in  the  density. 


Note 
The  system  of  integral  equations  (3)  was  given  independently  by 
BogolubovL-1  and  by  Born  and  Green'-  K     Closely  related  integral  equations 
had  previously  been  considered  by  Yvon  ^  -•  and  by  Kirkwood  and  Boggs  *■   ■'  • 
Mayer'-  •'  also  makes  extensive  use  of  integral  relations  between  the  dis- 
tribution functions. 
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